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ABSTRACT 


It has long been recognized that galactic cosmic rays are of such high energy that they tend to pass 
through available shielding materials resulting in exposure of astronauts and equipment within space 
vehicles and habitats. Any protection provided by shielding materials result not so much from stopping 
such particles but by changing their physical character in interaction with shielding material nuclei 
forming, hopefully, less dangerous species. Clearly, the fidelity of the nuclear cross-sections is essential 
to correct specification of shield design and sensitivity to cross-section error is important in guiding 
experimental validation of cross-section models and database. We examine the Boltzmann transport 
equation which is used to calculate dose equivalent during solar minimum, with units (cSv/yr), 
associated with various depths of shielding materials. The dose equivalent is a weighted sum of 
contributions from neutrons, protons, light ions, medium ions and heavy ions. We investigate the 
sensitivity of dose equivalent calculations due to errors in nuclear fragmentation cross-sections. We do 
this error ancdysis for all possible projectile-fragment combinations (14,365 such combinations ) to 
estimate the sensitivity of the shielding ccdculations to errors in the nuclear fragmentation cross-sections. 
Numerical differentiation with respect to the cross-sections will be evaluated in a broad class of materials 
including polyethylene, aluminum and copper. We will identify the most important cross-sections for 
further experimented study and evaluate their impact on propagated errors in shielding estimates. 

Introduction 

Particle transport equations are derived from continuum mechanics principles, (Wilson et al 1991). 
The particle flux in a shielding material is determined by balancing the change in particle flux across a 
small volume element of material with gains and losses caused by atomic and nuclear collisions within 
the material. The resulting equation is the well known Boltzmann equation in the continuous slowing 
down approximation given as 

[Q. • V - -j- ^ S j (E) + a j (E)]<pj (x. Q, E) = ^ J J dE ' dQ! a jk (E, E , Q., Q.' )<p k (x, Q.' , E ) ( 1 ) 

j k EQ. 

where </>j(x,Q.,E) represents the flux of type j-particles with atomic mass Aj at position x with motion in 
the direction Q. having energy E. Occurring in the equation (1) are the terms <7 j (E) which represents 

the macroscopic cross section, Sj(E) representing the linear energy transfer or change in energy per unit 
distance. The fragmentation of the projectile and target nuclei is represented by the fragmentation cross 
sections a jk (E,E',Q.,Q.') which represents the production cross section for type j particles with energy E 
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and direction El having a collision with a type k particle of energy E’ with direction El' . That is, <7 jk is a 

cross-section for producing ion j from a collision by ion k. These cross sections are composed of three 
parts and can be written 

a jk (E,E',£l.,£l') = <r k (E')v jk (E') f jk (E,E',El,El') (2) 

where V jk ( E' ) represents the average number of type j particles produced by a collision with a type k 

particle of energy E’. The term f jk (E,E',£l,Q.') is the probability density distribution for producing 

particles of type j of energy E into the direction El from the collision of the type k particle with energy E ’ 
moving in the direction El' . See Wilson, et al 1991. 

The propagation of galactic cosmic rays into shielding material (Cucinotta, et al 2003) is described by the 
above Boltzmann equation. The galactic cosmic ray (GCR) transport radiation code HZETRN, (Shin et al 
1992), was developed by NASA Langley Research Center around the 1980-1990 period. This code 
solves a one-dimensional form of the equation (1) and uses quality factors (ICRP, 1991) to calculate total 
radiation dose equivalent, in units of cSv/yr at various depths in a shielding material. The results are 
particularly important in determining radiation exposure of astronauts and electrical equipment on space 
missions. It has long been recognized that galactic cosmic rays are of such high energy that they tend to 
pass through available shielding materials resulting in exposure of astronauts and equipment within space 
vehicles and habitats. Any protection provided by shielding materials results not so much from stopping 
such particles but by changing their physical character by interaction with shielding material nuclei to 
hopefully less dangerous species. An understanding of these processes can be discerned by conducting 
various shielding simulations using the HZETRN computer code. Clearly, the fidelity of the nuclear 
cross sections are essential to correct specification of shield design and sensitivity of computational 
results to cross section error is important in guiding experimental validation of cross section models and 
the construction of cross section databases. 

In this paper we examine a one-dimensional form of the Boltzmann transport equation (1), which is 
solved by way of the HZETRN code subject to the galactic cosmic rays and the 1977 solar minimum 
environment. The HZETRN code calculates a radiation dose equivalent, with units of (cSv/yr), associated 
with various depths within a shielding material. The radiation dose equivalent is a weighted sum of 
contributions from protons, alpha particles, light ions, medium ions, heavy ions and neutrons using the 
ICRP-60 LET dependent quality factor (ICRP, 1991). We consider production terms arising from 170 
possible projectile-fragment cross-sections (Cucinotta, et al. 2003). The isotopes considered in the 
interactions are listed in the table 1 . The accuracy of the resulting dose equivalent is dependent upon the 
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projectile-fragment production cross-sections. In the Boltzmann equation the stopping power and 
nuclear total cross sections are known to a high degree of accuracy but the fragmentation cross sections 
are more dependent upon details of the nuclear processes and consequently exhibit more uncertainty when 
tested experimentally (Wilson, et al. 1995, Golovchenko et al. 2002). We investigate the sensitivity of 
the dose equivalent calculations due to production cross section errors using sensitivity analysis. We do 
this error analysis for all possible projectile-fragment combinations (14,365) to estimate the sensitivity of 
the shielding calculations to errors introduced into a single nuclear fragmentation cross-section. 
Numerical differentiation with respect to the cross sections is evaluated for the selected shield materials of 
polyethylene, aluminum and copper. We will identify the most important cross sections affecting the 
dose equivalent calculations. The purpose of the investigation is to evaluate the cross sections impact on 
propagated errors in shielding estimates associated with radiation protection studies. 

Error Analysis 

We introduce a small error into a single projectile-fragment cross section on the right-hand side of 
equation (1) and calculate the error in the calculated dose equivalent. Note that only one selected cross 
section out of 14,365 projectile-fragment cross sections on the right-hand side of equation (1) is being 
changed to calculate a resultant error in the dose equivalent. The error in the dose equivalent is obtained 
from a truncated Taylor series expansion in the parameter space defined by the set of 14,365 cross- 
section parameters. For example, in the case Ax is small one can use the first order error approximation 
for the error as 

~\r 

error = /(x + Ax, y,z,...)— /(x, y,z,...) = — Ax . 

3x 

The partial derivative term is taken as defining the sensitivity of the function f(x,y,z,...) in relation to 
changes in the single variable x. Here only a single cross section a ^ on the right-hand side of equation 
(1) is perturbed at a time to maintain parameter independence. These parameters appear on the right-hand 
side of equation (1) which is the gain side of the Boltzmann equation. The total cross section <jj and the 
stopping power Sj(E) on the left-hand side of equation (1), representing loss terms, has not been changed 
since these values are well known. One should note that this would technically violate the conservation 
of flux in the Boltzmann equation. This violation already exists resulting from uncertainty of especially 
the fragmentation cross sections as known from experiments (for example, Golovchenko et al. 2002). 

This paper evaluates the gradient of dose equivalent with respect to the fragmentation cross section 
parameter space at the location of nominal cross section values so that the effects of off-nominal values 
can be further investigated. Adding terms to the left-hand side of equation (1) to maintain a conservation 
law would constrain the off-nominal domain to selected trajectories in parameter space where flux is 
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conserved. However, the relative importance of the dose equivalent errors, obtained by this research, 
would remain the same (especially within the parameter subspace where flux is conserved). 

Consequently, we have made the assumption that by changing consecutively only one projectile- 
fragment cross-section at a time, out of 14,365 cross sections occurring on the right-hand side of equation 
(1), the corresponding error in the total dose equivalent is small (crudely, 1 part in 14,365) and the cross 
term variations are of even higher order and negligible ( similarly, 2 parts in 14,365). It should be further 
noted that the dose equivalent is a function of both energy and square of the charge of the ionizing 
radiation. It is proportional to the square of the charge of the nucleus. In the example given by 
56 Fe — > 55 Mn+p, the dose deposited before the fragmentation occurs is proportional to 56 2 =3 136, and 
after the fragmentation it is proportional to 55 -t- 1 =3026 . So the energy deposit is 110/3026 times larger 
before the fragmentation than after. It is known that such fragmentation processes are used to reduce the 
radiation damage effects. In this example, we have not taken into account the energy effects because for 
the high energy processes that this paper discusses the energies of the initial ionizing ion are very 
comparable to the fragmenting nuclei. 

We calculate the above Taylor series representation of the errors for each of all possible projectile- 
fragment combinations (14,365). These individual errors are then used to estimate the sensitivity of the 
shielding calculations to errors in all the nuclear fragmentation cross-sections. Numerical differentiation 
with respect to the cross sections is evaluated for the selected shield materials of polyethylene, 
aluminum and copper. We will identify the most important cross sections and corresponding ions for 
further experimental study. The purpose of the investigation is to evaluate the cross sections impact on 
propagated errors in shielding estimates associated with radiation protection studies. 

Unitarity Considerations 

We examine the role of unitarity in the evaluation of the sensitivity coefficients. We note that given a 
cross section data set [a jk + e jk ) where e jk represents the uncertainty in cross section knowledge, there is 

no a priori prescription for satisfying unitarity and one uses arbitrary models. If we are to arbitrarily 
apply a unitarity requirement then the appropriate functionality is 

H(v ik +e jk ) >H\U(a jk + £ j k )] = H{c jk +£ jk ) + Y J H((T, k + %') (3) 

where e Ik ' represents the modifications due to the chosen unitarity process denoted by the operator U and 
the sum extends over all l,k modified by the unitarity requirement. We will examine the relationship of 
variations of H(c t jk +e jk ) relative to variations in H[U(o jk + e jk )] . This relationship is clarified by using 

the chain rule of partial differentiation. Hence, one can write 
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D H[U(c T jk +£j k y\ = 


dH(a jk + e jk ) ^ y, dH(a lk +e, k ’) de lk ' 
tejk ^ tejk 


where D denotes a total derivative operator (both explicit and implicit differentials) and the sum extends 


over all parameters £ jk demanded by U. One form of unitarity is that any perturbation of a 


fragmentation cross section needs to be balanced by a change in the projectile absorption cross section o k 
as given by 

Z Aj (Oj k + £ jk ) = A k (Ok + e' k ) (5) 

where the sum extends over all the j fragments and only one value of £j k is nonzero in that description. In 
this case, e'k is numerically Aj £jk/Ak. An alternate unitarity principle by Dr. J. Letaw from the literature, 

( Letaw 1985), is to assume the balance is made in the lightest fragment as 

Ai (oik + £' ik) + Z Aj (Ojk + £ jk ) = A k o k (6) 

where the sum extends over j > 1 and only one value of £j k is specified as nonzero but for arbitrary j 
greater than one. In this case, £'ik is numerically -Aj £jk/Ai. It is clear that any arbitrary form of 
normalization for mass conservation can be used resulting in a broad range of arbitrary results for the 
“constrained sensitivity coefficients.” In any case, equation (5) is a poor choice as the absorption cross 
sections are accurately known for which the implied uncertainty factor £' k is misapplied. For a third 
alternate choice of renormalization of fragmentation models, that was used to improve comparisons with 
Ne beam experiments to conserve mass and charge, was applied equally to all fragments, see Wilson et 
al. (1991). We will now examine the impact of the two schemes of renormalization given by equations 
(5) and (6) on the sensitivity coefficients outcomes. 

Aside from factors dependent on kinematic variables, depth, and intensity at the boundary, the 
contribution to dose equivalent resulting from a given k primary ion is approximately proportional to 

H k = ZQj Zj 2 o jk + Q k Z k 2 ( 1 - o k ) (7) 

where the sum is over all j and Qj is a representative value of the quality factor. The total H is the sum 
over all j weighted by the k fluence at the boundary F k . It is clear from equation (7) that the e jk sensitivity 

(explicit derivative of H), apart from weighting with the k ion fluence at the boundary, is 

|^-=Q,Z , 2 (8) 

where we note that this is a simple explicit partial derivative as appears in the Taylor expansion above. 
The sensitivity as argued by renormalization requirements needs added contributions from implicit 
differentials related to the constraints. Renormalization suggests requiring one to evaluate the total 
derivative consisting of explicit and implicit variations as 
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( 9 ) 


, = Qj Zj 2 - (Aj/Ak) — = Qj Zj 2 - (Aj/Ak) Q k Z k 2 

total o<7 j 

where Aj /Ak is the derivative of e'k with respect to £jk as required by equation (5). Note that the main 
contributions from j ion fragments will be for Zk ~ Zj for which Ak ~ Aj and the total derivative is near 
zero implying little sensitivity under this renormalization rule. There is no great mystery as to why this 
occurs. If one always balances the mass loss in a given ion by an ion of nearly equal charge and mass, the 
resulting change in H is quite small. Only those cross sections for which Zj « Zk and Aj « Ak will 
appear as important as shown by the analysis in equation (9) using constraint equation (5). Applying the 
same process to the Letaw renormalization of equation (6) one attains 

Pj ZJ 

, = Qj Zj 2 - (Aj/AO- = Qj Zj 2 - (A j/Aj) Qj Zj 2 . (10) 

total d<T ljt 

For the galactic cosmic ray HZE ions, Aj ~ 2 Zj and the Qk is on the order of 20 to 30 while Qi is unity. 
It is clear that the derivative is large and nearly equal to the model independent sensitivity coefficient 
when Zk » Zi and near zero as Zk => Zj. What is clear from this simple analysis is that enforced 
different unitarity requirements introduces broad and conflicting results for the “sensitivity” as seen in 
the two examples above. The different choice of renormalization such as Letaw ‘s unitarity in equation 
(6) completely contradicts the choice of equation (5). Since the “unitarity requirement” is arbitrary and 
dependent on the choice of the unitarity renormalization scheme, we mainly focus on the unconstrained 
derivatives as best indicators of effects of cross section improvements. In all cases it is clear that the 
fundamental role is played by the sensitivity coefficients given by the explicit derivatives, but application 
of an arbitrary unitarity principle can lead to confusion in a final interpretation. 

Unitarity is a constraint that can be implemented in an infinite number of different ways some of 
which can highly bias the results of a sensitivity analysis as previously discussed. We have compared the 
results obtained in this study, where unitarity is not fully enforced, with selected numerical experiments 
where a form of unitarity is imposed. The unitarity constraint imposed is based upon the fact that total 
absorption cross sections are known with much greater accuracy than partial fragmentation cross sections. 
Hence we selected a normalization constraint that keeps the total absorption cross section constant while 
only changing the fragmentation cross sections. Unitarity can then be enforced by introducing an error in 
one fragmentation cross section a JK while simultaneously lowering the cross sections a Jk , for 
k=J+l,. . .,170 and k ^ K. We lower these fragmentation cross sections by a constant factor in order that 

170 

the sum ^<7 jk maintains its original value. We think that this method is less likely to introduce bias into 

*= 7+1 
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the calculation because the change in cross section is spread out equally throughout all the different 
fragmentation channels (similar to the Ne beam renormalization comparisons as shown by Wilson et al. 
1991) rather than arbitrarily choosing a few channels to change. We then performed numerical 
experiments where the dose equivalent was calculated for errors in J,K values associated with the top 5 
maximum errors listed in the table 3. The results from these selected numerical experiments was that the 
dose equivalent errors calculated with the imposed unitarity had the same values as previously obtained 
using our original calculations where renormalization was not used. 

Taylor series approximation of Errors 

Let H denote the total dose equivalent calculated by the HZETRN code using nominal cross sections 
from NUCFRG2 (Wilson et al. 1995). If a perturbation error £ t] > 0 is introduced into a single 

fragmentation cross sections rx ;; used to calculate of the total dose equivalent, then what kind of errors 

can be expected in the predicted radiation dose equivalent? By employing a Taylor series expansion one 
can show that 

3 // 

H( < 7 , + £ u )-H( <j, )=YV £.. + higher order terms. (11) 

ij^a 

Note that the matrix of partial derivatives on the right-hand side of equation (11) expresses the usual 
sensitivity of the dose equivalent H on the fragmentation cross section parameters . We neglect the 

higher order terms of the Taylor series expansion and consider a single error associated with an i,j 
combination where i and j have fixed values. If the first order term is large, then the higher order terms 
would need to be examined. By introducing an error into a single cross section we can obtain a numerical 
approximation for just one term from the right-hand side of equation (5). By going through all i, j 
combinations one obtains a numerical approximation for each term on the right-hand side of the Taylor 
series. We can then use these numerical approximations for the derivatives to formulate a Monte Carlo 
simulation given by 

H{Gij + Sijij y) = H(a ij ) + (12) 

i j <’ 

where //, ; are random variates from a normal distribution In this way one can model the effect of errors in 
calculating the total dose equivalent. 

The figure 1 illustrates the resulting Taylor series approximation for the errors in the dose equivalent 
at various depths within an aluminum shield due to the introduction of a single error in a specific k- 
projectile, j-fragment cross section without using any unitary considerations. For each given j,k value the 
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cross section <7 jk was replaced by 2 <J jk in order to control effects of round-off error. This was done for 

all 14,365 combinations of j and k. Also illustrated in the lower panel of these figures is the summation 
of fragment errors associated with a given projectile k. The figures 2 and 3 are similar figures for copper 
and polyethylene respectively. 

Results and Discussion 

Typical dose equivalent constituents are illustrated in the figures 5 and 6. The figures 1,2, and 3 
illustrate the errors obtained in the dose equivalent calculations. These errors are a numerical 
representation of the derivative times the error representing a single term from the right-hand side of 
equation (9). The highest five peaks occurring in the figures 1,2 and 3, ordered from highest to lowest, 
correspond to the isotopes of 2 6Fe , 8 0 , wSi , i 2 Mg and 6 C . These elements are the five most 
abundant heavier elements in the Galactic cosmic ray composition. The dose equivalent errors for 
aluminum, copper and polyethylene are more pronounced along the line of lower J-values and along the 
line where j ~ k . This is due to projectile interaction with shield material which causes few nucleon 
removal followed by resultant particle movement into the material showing up as the errors near the 
j = k line. Consequently, the dose equivalent errors correspond to few nucleon removal followed by 
projectile fragment continuation after interaction. 

Monte Carlo studies which incorporate errors from all sources were conducted by generating random 
values for tj from a normal distribution with mean // = 0 , and standard deviation, std=0.75. We then 

conducted 50 Monte Carlo simulations for errors produced in the dose equivalent at various depths. 
Results from these simulations indicate low overall errors produced in the dose equivalent calculations. 
These results are indicated by the error bars on the curves given in figure 4. Also indicated in the figure 4 
are curves obtained where all cross sections were set equal to zero. The tables 2,3 and 4 give the fifty 
maximum errors associated with 20 g/cnT of shield material for the materials of aluminum, polyethylene 
and copper respectively. 

Similar sensitivity studies for an aluminum shield can be found by Townsend et al 1992, where in that 
study the production cross-sections <J jk where replaced by per k . where p varied from 0.5 to 1.5 and the 

dose equivalent calculations where compared with the nominal values obtained when p=l. 

One might conclude that the results of this study suggest that for the 1977 solar minimum environment 
and GCR environment, the dose equivalent is not very sensitive to errors in the production cross sections. 
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However, one can not conclude that secondary particle production is not important in general. Note that 
the production cross section errors can greatly affect the shielding requirements for a given dose limit. 
This is illustrated by the horizontal lines depicted in the figure 4. Particle production from lighter nucleon 
removal has the largest effect on dose equivalent from the HZETRN code. The largest change in dose 
equivalent is due to light, medium and heavy ion drop off with depth into shielding material. Emphasis 
on the most abundant elements in the Galactic cosmic rays and their particle interaction with the shield 
material is the most important issue in determining dose equivalent predictions. Future efforts are being 
directed toward developing additional physics and mathematics for describing the propagation of particles 
through a shielding medium. 

Acknowledgement: This research was sponsored by NASA grant NAG- 1-03075. 
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Figure 2. Error and cumulative error graphs for 
Copper shield at various depths 
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Figure 3. Error and cumulative error graphs for 
polyethylene shield at various depths 
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Figure 4. Dose equivalent comparison with and without 
production cross-sections equal to zero. 
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Figure 6. 

Dose contribution polyethyelene 
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Range 

J 

K 

Error 

Fragment 

Projectile 


Range 

J 

K 

Error 

Fragment 

Projectile 

1 

162 

164 

2.52E-01 

2 6 Fe 55 

Fp 56 
26^ e 


26 

123 

164 

5.55E-02 

2 oCa 43 

Fp 56 

26™ 

2 

161 

164 

2.09E-01 

25 M 11 

Fp 56 
26^ e 


27 

133 

164 

5.46E-02 

22 Ti 46 

Fp 56 
26^ e 

3 

159 

164 

1.89E-01 

25 M 11 54 

Fp 56 
26 ^ e 


28 

126 

164 

5.39E-02 

2 ,Sc 44 

2 6 Fe 56 

4 

156 

164 

1.33E-01 

25 M 11 

2 6 Fe 56 


29 

21 

23 

5.31E-02 

6C“ 

6 c 12 

5 

33 

35 

1.29E-01 

8 o 15 

8 o 16 


30 

49 

58 

5.09E-02 

ioNe 21 

i 2 Mg 24 

6 

153 

164 

1.21E-01 

2 4 Cr 52 

Fp 56 

26 ™ 


31 

140 

164 

5.06E-02 

V 48 
23 V 

Fp 56 
26^ e 

7 

150 

164 

1.12E-01 

24 Cr 51 

2 6 Fe 56 


32 

120 

164 

5.02E-02 

2 oCa 42 

Fp 56 
26^ e 

8 

2 

35 

1.11E-01 

,H> 

8 o 16 


33 

52 

58 

4.80E-02 

uNa 22 

12 Mg 24 

9 

68 

71 

1.07E-01 

i 4 Si 27 

i 4 Si 28 


34 

147 

164 

4.59E-02 

rV° 

24^r 

Fp 56 

26™ 

10 

143 

164 

9.86E-02 

23' V 49 

Fp 56 

26™ 


35 

113 

164 

4.45E-02 

T^ 40 

1 9 K 

Fp 56 

26™ 

11 

55 

58 

9.57E-02 

i 2 Mg 23 

, 2 Mg 24 


36 

160 

164 

4.37E-02 

2 6 Fe 54 

2 6 Fe 56 

12 

2 

23 

9.51E-02 

1H 1 

6 c 12 


37 

26 

35 

4.35E-02 

6 c 13 

8 o 16 

13 

67 

71 

9.18E-02 

13AI 27 

i 4 Si 28 


38 

44 

47 

4.34E-02 

ioNe 19 

ioNe 20 

14 

29 

35 

8.83E-02 

7 N 14 

8 o 16 


39 

154 

164 

4.22E-02 

25 M 11 52 

Fp 56 
26^ e 

15 

136 

164 

8.65E-02 

T*47 

22 A 1 

2 6 Fe 56 


40 

116 

164 

4.16E-02 

9 o 18 

Fp 56 
26^ e 

16 

146 

164 

8.58E-02 

23 V 50 

2 6 Fe 56 


41 

132 

164 

4.04E-02 

2 ,Sc 46 

2 6 Fe 56 

17 

32 

35 

8.17E-02 

7 n 15 

8 o 16 


42 

23 

35 

3.85E-02 

6 C 12 

8 o 16 

18 

54 

58 

7.60E-02 

uNa 23 

, 2 Mg 24 


43 

105 

164 

3.81E-02 

isAr 38 

2 6 Fe 56 

19 

129 

164 

7.05E-02 

21 SC 45 

2 6 Fe 56 


44 

54 

71 

3.77E-02 

uNa 23 

i 4 Si 28 

20 

1 

35 

6.97E-02 

n 

8 o 16 


45 

64 

71 

3.65E-02 

13AI 26 

i 4 Si 28 

21 

60 

71 

6.28E-02 

12 Mg 25 

i 4 Si 28 


46 

17 

23 

3.64E-02 

5 B 10 

6 c 12 

22 

139 

164 

6 . 11 E -02 

rp-48 

22 A 1 

2 6 Fe 56 


47 
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153 

3.58E-02 

24 Cr 5 1 

24 Cr 52 

23 
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5.97E-02 
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6 c 12 
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111 

114 

3.44E-02 

2 oCa 39 

Pa 40 

20 '-& 

24 

155 

164 

5.90E-02 

24 Cr 53 

2 6 Fe 56 


49 

136 

139 

3.40E-02 

2 2 Ti 47 

Ti 47 
22 A 1 

25 

63 

71 

5.85E-02 

12 Mg 26 

i 4 Si 28 


50 

51 

58 

3.39E-02 

ioNe 22 

12 Mg 24 


Table 2. Dose equivalent sensitivity with top 50 cross sections 

Maximum errors (cSv/yr) associated with 20 g/cm 2 of aluminum shield. 


15 



Range 

J 

K 

Error 

Fragment 

Projectile 


Range 

J 

K 

Error 

Fragment 

Projectile 

1 

162 

164 

2.51 E-01 

26 Fe 55 

26 Fe 56 


26 

157 

164 

6.65E-02 

26Fe 

26Fe 

2 

33 

35 

2.06E-01 

8 o 15 

8 o 16 


27 

155 

164 

6.59E-02 

r r 53 

24^4 

c _56 

26Fe 
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159 

164 

1.82E-01 

25Mn 

1 -^ 56 

26Fe 


28 

133 

164 

6.49E-02 

Ti 46 
22 1 1 

C -56 

26Fe 

4 

161 

164 

1.76E-01 

25Mn 

26Fe 


29 

147 

164 

6.01 E-02 

24 Cr 50 

26 Fe 56 

5 

153 

164 

1.68E-01 

24 Cr 52 

26 Fe 56 


30 

17 

23 

5.78E-02 

5 B 10 

6 c 12 

6 

29 

35 

1.54E-01 

7 N 14 

8 0 16 


31 

20 

23 

5.70E-02 

5 B 11 

6 c 12 

7 

55 

58 

1.32E-01 

i 2 Mg 23 

i 2 Mg 24 


32 

26 
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/-\13 
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68 
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1.30E-01 

i 4 Si 27 

i 4 Si 28 


33 

64 
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i 3 ai 26 

i 4 Si 28 
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32 

35 

1.27E-01 

7 N 15 

8 0 16 


34 

120 

164 

5.43E-02 

20 Ca 42 

26 Fe 56 

10 

150 

164 

1.27E-01 

2 4 Cr 51 

26Fe 


35 

49 

58 

5.39E-02 

ioNe 21 

i 2 Mg 24 

11 

156 

164 

1.14E-01 

25Mn 

c _56 

26Fe 


36 

63 

71 

5.33E-02 

i 2 Mg 26 

i 4 Si 28 

12 

143 

164 

1.1 IE-01 

23 V 49 

26 Fe 56 


37 

123 

164 

5.23E-02 

20 Ca 43 

26 Fe 56 

13 

21 

23 

1.05E-01 

eC 11 

6 c 12 


38 

44 

47 

5.17E-02 

ioNe 19 

ioNe 20 

14 

2 

35 

9.97E-02 

1H 1 

8 0 16 


39 

1 

35 

5.02E-02 

n 

8 0 16 

15 

136 

164 

9.92E-02 

22 Ti 47 

26 Fe 56 


40 

140 

164 

4.91 E-02 

23 V 48 

r-^56 

26Fe 

16 

67 

71 

9.1 3E-02 

13AI 27 

i 4 Si 28 


41 

126 

164 

4.85E-02 

21 Sc 44 

26 Fe 56 

17 

160 

164 

9.03E-02 

26 Fe 54 

26 Fe 56 


42 

47 

58 

4.58E-02 

ioNe 20 

i 2 Mg 24 

18 

2 

23 

8.67E-02 

1 H 1 

6 c 12 


43 

154 

164 

4.57E-02 

25 Mn 52 

26 Fe 56 

19 

54 

58 

8.39E-02 

nNa 23 

i 2 Mg 24 


44 

58 

71 

4.55E-02 

i 2 Mg 24 

i 4 Si 28 

20 

23 

35 

8.37E-02 

6 C 12 

8 0 16 


45 

150 

153 

4.32E-02 

r r 51 
24'-' • 

24 Cr 52 

21 

139 

164 

7.76E-02 

22 Ti 48 

i -^5 6 

26Fe 


46 

1 

23 

4.30E-02 

n 

6 c 12 

22 

129 

164 

7.75E-02 

21 Sc 45 

56 

26Fe 


47 

23 

29 

4.13E-02 

6 c 12 

7 N 14 

23 

146 

164 

7.17E-02 

23 V 50 

26 Fe 56 


48 

43 

47 

3.91 E-02 

9 F 19 

ioNe 20 

24 

52 

58 

7.05E-02 

nNa 22 

i 2 Mg 24 


49 

136 

139 

3.88E-02 

Ti 47 
22 1 1 

Ti 48 
22 1 1 

25 

60 

71 

6.86E-02 

i 2 Mg 25 

i 4 Si 28 


50 

81 

84 

3.81 E-02 

C 31 

16 ^> 

r »32 

16 ^> 


Table 3. Dose equivalent sensitivity with top 50 cross sections 

Maximum errors (cSv/yr) associated with 20 g/cm 2 of polyethylene shield. 
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Rank 

J 

K 

Error 

Fragment 

Projectile 


Rank 

J 

K 

Error 

Fragment 

Projectile 

1 

162 

164 

2.90E-01 

2 6 Fe 55 

Fp 56 

26 ^' e 


26 

133 

164 

4.09E-02 

22 Ti 46 

Fp 56 

26 ™ 

2 

161 

164 

1.87E-01 

25M11 

Fp 56 

26 ™ 


27 

126 

164 

4.08E-02 

2 ,Sc 44 

26 Fe 56 

3 

159 

164 

1.32E-01 

25 MlV 4 

Fp 56 

26 ™ 


28 

63 

71 

3.90E-02 

, 2 Mg 26 

i 4 Si 28 

4 

156 

164 

9.42E-02 

25M1I 

2 6 Fe 56 


29 

120 

164 

3.86E-02 

Tr, 42 

20*- a 

Fp 56 

26 ™ 

5 

2 

35 

8.84E-02 

1H 1 

sO 16 


30 

140 

164 

3.75E-02 

23 V 48 

26 Fe 56 

6 

153 

164 

8.67E-02 

2 4 Cr 52 

2 6 Fe 56 


31 

150 

153 

3.73E-02 

24 Cr 5 1 

24 Cr 52 

7 

33 

35 

8.51 E-02 

uNa 22 

8 o 16 


32 

113 

164 

3.45E-02 

1^40 

1 9 K 

Fp 56 

26 ™ 

8 

150 

164 

8.1 IE-02 

24 Cr 51 

2 6 Fe 56 


33 

136 

139 

3.37E-02 

22 Ti 47 

22 Ti 48 

9 

68 

71 

7.86E-02 

i 4 Si 27 

i 4 Si 28 


34 

21 

23 

3.37E-02 

6C 11 

6 c 12 

10 

2 

23 

7.66E-02 

,H‘ 

6 c 12 


35 

49 

58 

3.34E-02 

ioNe 21 

12 Mg 24 

11 

67 

71 

7.24E-02 

13AI 27 

i 4 Si 28 


36 

147 

164 

3.34E-02 

f> 50 

24 ^r 

Fp 56 

26 ™ 

12 

143 

164 

7.24E-02 

23 V 49 

Fp 56 

26 ^' e 


37 

116 

164 

3.19E-02 

ir 4 1 

i 9 k 

Fp 56 

26 ™ 

13 

55 

58 

6.70E-02 

oMg 23 

, 2 Mg 24 


38 

111 

114 

3.19E-02 

20 Ca 39 

fr, 40 

20*-'* 

14 

136 

164 

6.45E-02 

T '47 

22 A 1 

Fp 56 

26 ™ 


39 

52 

58 

3.14E-02 

..Na 22 

12 Mg 24 

15 

146 

164 

6.22E-02 

V 50 
23 V 

2 6 Fe 56 


40 

160 

164 

3.07E-02 

26 Fe 54 

Fp 56 

26 ™ 

16 

1 

35 

5.76E-02 

n 

8 o 16 


41 

159 

161 

3.05E-02 

2 5Mn 54 

hMii 5 

17 

29 

35 

5.65E-02 

7 n 14 

8 o 16 


42 

154 

164 

3.04E-02 

25 Mn 52 

Fp 56 

26 ™ 

18 

54 

58 

5.65E-02 

uNa 23 

, 2 Mg 24 


43 

132 

164 

3.01 E-02 

21 SC 46 

Fp 56 

26 ™ 

19 

32 

35 

5.38E-02 

7 N 15 

8 o 16 


44 

105 

164 

2.97E-02 

isAr 38 

26 Fe 56 

20 

129 

164 

5.33E-02 

21SC 45 

2 6 Fe 56 


45 

44 

47 

2.85E-02 

ioNe 19 

ioNe 20 

21 

1 

23 

5.05E-02 

n 

6 c 12 


46 

2 

58 

2.83E-02 

1H 1 

12 Mg 24 

22 

139 

164 

4.50E-02 

rp -48 

22 A 1 

2 6 Fe 56 


47 

26 

35 

2.78E-02 

6 c 13 

8 o 16 

23 

60 

71 

4.24E-02 

, 2 Mg 25 

i 4 Si 28 


48 

149 

153 

2.70E-02 

V 51 

23 V 

24 Cr 52 

24 

123 

164 

4.21 E-02 

2 0 Ca 43 

Fp 56 

26^' e 


49 

81 

84 

2.62E-02 

i 6 s 31 

e 32 

16 ^ 

25 

155 

164 

4.17E-02 

24 Cr 53 

Fp 56 

26 ™ 


50 

54 

71 

2.57E-02 

..Na 23 

i 4 Si 28 


Table 4. Dose equivalent sensitivity with top 50 cross sections 

Maximum errors (cSv/yr) associated 20 g/cm 2 copper shield. 
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